Relaxation times hierarchy in two-component quasiparticle gas 
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sir.. 



A quasiparticle description of various condensed media is 
a very popular tool in study of their transport and thermo- 
dynamic properties. I present here a microscopic theory for 
the description of diffusion processes in two-component gas 
of quasiparticles with arbitrary dispersion law and statistics. 
Particularly, 1 analyze the role of interaction within each sub- 
system (i.e. between identical quasiparticles) in relaxation of 
the whole system. The approach for solution of such kinetic 
problem allows to study the most important limiting cases and 
to clarify their physical sense. Classical results for diffusion 
coefficient of light particles in a massive gas (Lorentz model) 
and of massive particles in a light gas (Rayleigh model) are 
obtained directly from the general solution without using ar- 
tificial approaches, as it was done earlier. This particularly 
provide a possibility to generalize these popular models on 
quasiparticle systems. 

63.20.-e, 66.10.Cb, 66.30.-h, 66.90.+r, 51.10.+y 



I. INTRODUCTION 

It is well known that variety of properties of some 
condensed media can be described by the interaction 
processes in quasiparticle gases. These are, say, trans- 
verse and longitudinal phonons in solids, phonons and 
magnons in magnetically ordered materials, phonons and 
rotons in supcrfruid helium, conduction electrons and 
holes in semiconductors etc. As a rule, to study the 
dissipative properties of such systems investigators use 
either classical kinetic theory in its simplest form or 
some semi-intuitive models, which lead sometimes to 
quite ambiguous results in case of quasiparticle systems. 
Here I present the theory of diffusion processes in two- 
component quasiparticle systems, which in general is in- 
dependent on particular quasiparticles statistics and dis- 
persion law. Let me first briefly review the present state 
of classical kinetic theory and analyze its pitfalls in de- 
scription of quantum quasiparticle systems of condensed 
media. 

The kinetic theory of gases in modern understanding 
takes its beginning from pioneering work by MaxwellEJ, in 
which he has proved the law for distribution of velocities 
of molecules in homogeneous equilibrium gas (so-called 
Maxwell velocities distribution) and law of an equidis- 
tribution of average energy of molecules in a mixture of 
gases. These his results were updated and improved in 
further works devoted to the theory of inhomogeneous 
gases (about a history of a problem see Ref. 2). How- 



ever as basis of all mathematical methods of the modern 
kinetic theory it is necessary to consider apparently the 
basic works by BoltzmannH, in which both the H-theorem 
was proved and the classical Boltzmann equation was in- 
troduced. 

Boltzmann equation is the integro-differential equation 
describing the collisional behaviour of rarefied gas. Un- 
til now it remains to be a basis of the kinetic theory of 
gases and it appears to be very fruitful not only for a re- 
search of classical gases, which Boltzmann himself kept in 
mind, but - with an appropriate generalization - for study 
of electron transport in solids and plasma, transport of 
neutrons in nuclear reactors, phonon and roton transport 
in superfluid liquids, and transport of a radiation in at- 
mospheres of stars and planets. For past 130 years these 
researches have led to significant achievements both in 
new areas, and in old one. 

Generally, kinetic equation of the Boltzmann type (the 
equation describing the evolution of a single-particle dis- 
tribution function in the phase space) represents the 
integro-differential equation, which remarkable property 
is the nonlincarity of collision term. Just this fact makes 
a main obstacle in construction of methods, for solution 
of the kinetic equation. The monographiesa'tl are devoted 
to a detail exposition of such methods in case of classical 
gas systems. 

In majority of experimental problems there is no ne- 
cessity to use the detail microscopic description of gas 
systems at a level of distribution functions. As a rule, 
investigation of physical processes in macroscopic sys- 
tems is carried out at a less detail level of hydrodynamic 
variables. Since these variables are determined through 
the moments of a distribution function, then, as a rule, 
the detail study of main moments of distribution func- 
tion appropriate to collisions invariants is required, but 
not the distribution function itself. Thus, the connection 
between the kinetic theory and hydrodynamics appears 
to be one of the main problem. In particular, one of 
the main aspects of this problem is the determination of 
transport coefficients, such as diffusion coefficient, vis- 
cosity (first and second), thermal conductivity appeared 
in equations of hydrodynamics of a viscous liquidO. 

In spite of so long history of physical kinetics, today we 
have rather small number of approaches to a solution of 
the kinetic equations. All these approaches and methods 
were formed depending on concrete problems, on which 
were directed. Most general methods of research of non- 
equilibrium state of classical and quantum gases were 
directed rather on demonstrating of a mathematical re- 
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solvability (or insolubility) of certain basic problems in 
principle, than on construction of the serviceable theo- 
ries, suitable for a solution of concrete physical problems. 

The classical methods of a solution of the kinetic equa- 
tions allow to derive the kinetic coefficients as series ex- 
pansions on an infinite set of orthogonal polynomials. 
However, it appeared to be very difficult to use these 
classical results for numerical calculation and analysis of 
physical processes in real systems. In particular, it is 
caused by impossibility of selection of the contribution 
from different types of interactions to various kinetic co- 
efficients, while many physical systems behave qualita- 
tively differently for different ratio between a speed of 
a relaxation inside each subsystem and between subsys- 
tems. 

In classical gases with point-like interaction this prob- 
lem is less important because the speed of relaxation in- 
side subsystems unambiguously determined by their mass 
and concentration. In the quantum case, when we talk 
about quasiparticles, the situation becomes more compli- 
cated. The mechanism of interaction between quasipar- 
ticles can be any independently on simple macroscopic 
parameters. For example, in some cases such notation as 
mass cannot be well defined at all (say, what is the mass 
of phonon?). 

The mathematical theory of transport processes is 
most sequentially advanced for mixtures of classical 
gases, the evolution of which is described by a set of 
Boltzmann equations. A basis of classical methods for 
solution of the Boltzmann equation in case of one 
component gas is the formal expansion of distribution 
function in power series of some parameter er in the form 
/ = /(°) + o-/« + o- 2 / (2) + • • so that function /(°) 
corresponds to statistical equilibrium. In this case pa- 
rameter a is some scale factor for density, the physical 
sense of which can be different depending on a particu- 
lar problem. As a rule, this parameter is formally con- 
sidered to be small, so that the solution of the kinetic 
equation represents a problem of singular perturbation! 
The most successful methods of a solution of the kinetic 
equation, such 3*-fhc Hilbert method and the Chapman- 
Enskog methodatl, are based on this principle. In spite 
of the successes of the Chapman-Enskog method in the 
description of connection of the kinetic theory with the 
equations of hydrodynamics (the Navier-Stokes equations 
appear already in the first order in parameter a), the 
explicit expressions for kinetic coefficients have a rather 
complicated form. Main defect of these expressions is 
that already in the first order in parameter a their anal- 
ysis becomes practically impossible. The situation be- 
comes more problematic in the case of two-component 
gas. The infinite series of integral brackets containing 
Sonin polynomials does not allow to select explicitly the 
contributions from different types of interactions in a sys- 
tem to various dissipative coefficients. This frequently 
leads to necessity to use various ungrounded approxima- 
tions, such as the Charoman-Cowling approximation^ or 
Kihara approximationB. 



At the same time, in spite of some successesB, there is 
no yet consistent mathematical theory for deriving the 
dissipative coefficients in gases of quasiparticles. And, 
naturally, the problem of distinguishability of contribu- 
tion from interactions between identical quasiparticles 
and between different subsystems has not been solved 
yet. This problem has his own history. In practice, 
when analyzing particular physical systems, many physi- 
cists use some model approximations for collision inte- 
grals. The roost popular one is the so-called BGK- 
approximatiorJj. In its simple form BGK-approximation 
leads sometimes to quite confusing results. Really, let us 
have several relaxation mechanisms in the system, and, 
therefore several characteristic times depending on mo- 
menta or energies of quasiparticles. The final observed 
quantities should be obtained by averaging this times 
in some manner. In particular problems the following 
question frequently appears: What must be averaged, 
e.i. the time or the rate (inverse time)? and how to 
obtain the real relaxation time, i.e. by summation of 
the times or by summation of the rates? This uncer- 
tainty led to many confusing situations. For example, 
more than for 15 years there were two different theories 
for the mobility of two-dimensional electseai gas local- 
ized over the free surface of liquid heliurr£3. The first 
theoryM assumes the mobility to depend on the aver- 
aged characteristic time of electron-ripplon interaction 
(ripplons are the quantized surface waves of liquid he- 
lium). This theory well describes the experimental data 
for small electron densitytl Another theoryO assumes 
the mobility to be determined by the inverse averaged 
rate of the same interaction— These theoretical results 
well fit another experimentslij with large electron den- 
sity. The problem of relationship between these two re- 
sults has naturally appeared. Similar confusion—was in 
the theory of dissipative processes in superfluidsE-a The 
analogous situation took placeJpr some time in theory of 
thermal conductivity in solidal3. 

The aim of this paper is to present alternative approach 
for solution of the system of linearized kinetic equations 
for two-component gas of quasiparticles with arbitrary 
statistics and dispersion law. The theory explicitly ac- 
counts the all types of interactions in the system. This 
allows to analyze the contribution of interaction between 
identical particles into the relaxation of the whole sys- 
tem. I do not restrict myself here in the frame of partic- 
ular system. So, the results obtained here can be applied 
to any quantum system, whose dissipative properties are 
determined by the processes in two-component gas of 
quasiparticles. 

The paper is constructed as follows. In Section II I 
formulate the problem mathematically and carry out the 
linearization procedure. Section III is devoted to the pro- 
cedure of inversing of collision operator by the use of pro- 
jection operators method. In Section IV I derive the ex- 
act solution for characteristic diffusion time and analyze 
all limiting cases. The classical systems of Lorentz and 
Rayleigh gases and their generalization to the quasipar- 
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ticle systems are considered in Section V. In Section VI I 
consider the generalization of the Kihara approximation 
on the quasiparticle quantum systems. The outlines and 
conclusions are given in Section VII. 



II. GENERAL EXPRESSIONS 

So, consider stationary nonequilibrium state of gaseous 
mixture of quasiparticles of two species. The most inter- 
esting relaxation process in such a system is a diffusionllZI, 
so I will concentrate on the diffusive processes. 

One of the most essential advantage of the offered the- 
ory undoubtedly is the fact, that it is correct for quantum 
gases with any statistics and any dispersion of quasipar- 
ticles. All main outcomes remain applicable for both sys- 
tems with nonzero chemical potential, and with chemical 
potential equal to zero (when the number of quasiparti- 
cles is not conserved). The constructed theory does not 
meet principal difficulties in generalization on multicom- 
ponent systems {i.e. on systems with number of active 
components exceeding two ones). 

The evolution of distribution functions fk (subscript 
k = 1, 2 numbers the component of the mixture) can be 
described by the following set of kinetic equations 



Vfc 7^ = °k]{fk 
3=1,2 



/,•; 



(fc = l,2), 



(1) 



where = de k / dp k is the velocity of quasiparticle of fcth 
type, efc and pt are its energy and momentum, respec- 
tively; r is the coordinate, Ckj{fk, fj) is the collision in- 
tegral, which is a functional of the distribution functions 
of mixture components 1 and 2. The particular form of 
these collision integrals depends on the concrete physical 
problem. To find the diffusion coefficient let us consider 
the stationary nonequilibrium state of two-components 
1 and 2 of the mixture, in which the quasiparticles num- 
ber densities are functions of coordinate r. In particular, 
for the gas of thermal excitations such situation can be 
realized by creation of constant temperature gradient. 

Under the considered conditions there are stationary 
gradients of partial pressure of components, which result 
in flow of quasiparticles. This flow is determined by mo- 
mentum current density 



Jfe 



3=1,2 F 



(* = 1,2), 



(2) 



where p k is the normal density of fcth component of the 
system, p = p\ + p2 is the total density, dkj is the matrix 
of diffusion times. The partial pressure of quasiparticles 
is determined in standard mannerO 



(3) 



where dTj is the measure in phase space. The density-i of 
fcth component can be written in the universal formlla 



(4) 



where f' k = dfk/de k . Relations (H)-(0) are suitable for 
the quasiparticles with arbitrary dispersion law, statis- 
tics, and chemical potential (I mean both for zero chem- 
ical potential and for nonzero one). Note, the definition 
of normal density (|J) does not depend explicitly on such 
notations as mass or number density. 

To exclude convective transport in a quasiparticle sys- 
tem, i.e. to investigate only dissipative processes, here- 
inafter I shall consider the sum of partial pressure of dif- 
ferent components of a mixture to be a constant, so put 
P = Pi + P-2 = const. The analytical relation between 
matrix elements of a matrix of diffusion times dkj and 
usual diffusion coefficient D of a binary mixture can be 
obtained by a direct comparison of expression (^|) with 
definition of diffusion coefficient. Thus, one can use the 
ordinary gas-dynamic definition of current density of fcth 
component of a mixture (see e.g. Ref. 5) 



VkfkdT 



k ■ 



(5) 



It appears that for various both quantum and classical 
physical systems diffusion, coefficient can be written in 
the most general formli3u3 



D = u d t d , 



(6) 



where «d is the characteristic velocity, whose analytic 
form depends on the particular dispersion law and statis- 
tics of quasiparticles, td is the characteristic diffusion 
time to be determined. 

According to the relations (^) and (|J) , for deriving the 
diffusion coefficient (Q) it is necessary to solve a system 
of the kinetic equations (1). Below for definiteness the 
diffusion in a system with conserved number of quasipar- 
ticles (and, therefore nonzero chemical potential) will be 
considered. The calculation for case with not conserved 
number of quasiparticles can be carried out within the 
similar framework. 

Since we are interested in the theory within linear re- 
sponse approximation, let us assume the deviation of dis- 
tribution functions fk from their locally equilibrium val- 
ues f^ to be small. So, put as usual 



fk 



AO) 

Jk 



5fk 



\Sfk\<f, 



(0) 



(7) 



The small deviation Sfk can be conveniently rewritten in 
the form 



Sfk = - 



-9k 



(8) 



with unknown functions g k . Linearizing the system of 
kinetic equations (|l|) we come to the system of linear 
integro-differential equations for unknown quantities g k , 
which determine the degree of system perturbation 
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Vfc dP k 
n k dr 



Ckk9k + C kj {g k + gj) 
(fe,j = l,2;fe ^ j). 



(9) 



Here n k is the number density of kth component, C k j are 
the linearized collision operators for the collisions within 
each component (k = j) and for the collisions between 
different quasiparticlcs (k ^= j). The acting of this opera- 
tors on an arbitrary function of momentum, say £(Pfc,i)i 
is determined by the particular form of collision integrals 
appeared in equation (|l|). If we deal with the ordinary 
binary collision integral with probability density function 
W k j(p k pj\p' k p'j) then we obtaint2l 

Cfcj£(Pfc,i) = J Wkj(PkPj\p'kPj)fj 0) (Pj) 

x{i±/i 0) (p fe )}"{i±/r(pD}{i±/i 0) (p;)} 

x [z(p' kd )-z(p k jj\dr J dr' k dr' j (10) 

for k 7^ j, and 

C kk t(Pk) = / Wkk(PkP\p' k P')fk \p) 

x{i±/r(p,)}"{i±/f(p',)}{i±/i o) (p')} 

x K( P ' fe ) + C(p') - t{Pk) - C(p)] drdr' fe dr' (11) 



for k = j . The plus and minus in eqs. ( |10| ) and (|llj) signs 
correspond to bosons and fermions, respectively. 

Note that in fact our general approach allows to ac- 
count in such a manner not only binary collisions, but 
a variety of more specific types of interaction, such as 
decay or conversion processeso, creation or annihilation 
of quasiparticles, interaction with boundaries and point 
defects etc. 



III. INVERTING OF COLLISION OPERATOR 

According to the relation the gradient of partial 
pressure for quasiparticles with nonzero chemical poten- 
tial can be written as 



dPk _ ( d^ k 
dr \ dr 



(12) 



T— const 



where \i k is the chemical potential of fcth subsystem, T 
is the mixture temperature. 

For further calculations it is convenient to present sys- 
tem (^J) in the compact matrix form 



fc=l,2 



(13) 



where 



IV'i) 



Vl?^ 





v 2 n 2 1 



1.9) = 



<Jl 

<J2 



(14) 



are the two-component ket- vectors, defined in infinite 
Hilbert space to be specified. The collisional operator 
matrix C can be conveniently written as a sum 



C = J + S, 
where the operator matrix 

t _ | C12 C12 
021 021 



(15) 



(16) 



contains only the quantities corresponding to interactions 
between quasiparticles from different subsystems, and 



S 



C u 
G 



■22 



(17) 



corresponds to the relaxation inside each subsystem. To 
define completely the Hilbert space I am working in, let 
us introduce the scalar product of two-dimensional bra- 
vector (C| = (Ci(Pi),C2(P2)| and ket-vector |x) = ((xl) f 
in the following manner 



(Clx) 



k=l,2 



(Cfelxfe) 



fe=l,2 



CkXk 



Oft 



(18) 



where (Q k \ and \xk) are the corresponding one 
component bra-vector and ket-vector, respectively. It is 
easy to verify that with such a choice of scalar product 
(U|) the collision operator C becomes hermitian. 

System ( |l3| ) is the system of non-uniform linear inte- 
gral equations. According to the general theoryof inte- 
gral equations the sought solution \g) of system (|lj) must 
be orthogonal to the solution of corresponding uniform 
equations 



C\<h 



0. 



(19) 



The normalized solution of equation ( |19[ ) can be written 
in the following form 



|0i 



1 



(20) 



This vector \cj>i) corresponds to the total momentum of 
our two-component quasiparticle system. In general, 
equation (|l^) has another solutions, corresponding to 
conservation of energy, particle number etc. I account 
here only the momentum conservation law because the 
sought solution can depend only on the quasiparticle mo- 
menta (see initial equations (|l3|)). It is convenient to 
write the formal solution of equation ( |l3| ) so that the 
orthogonality condition 
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0#i 



o 



(21) 



is contained explicitly in the solution. For this purpose 
let us define the projection operator V n onto the subspace 
orthogonal to the vector \<j>\), 



1 



\<h){<h 



(22) 



As a result, the formal solution of equation ([13]) can be 
written in the form 



dP k 



k=l,2 



(23) 



Further we must insert the solution (23) into the expres- 
sion for current density (|5|) keeping in mind the relations 
(0) and (|J). Comparing the obtained result with defini- 
tion (|^) for a matrix of diffusion times we come to 



P2 



TD, di2 = G?21 



TD, 



t22 



Pi 



Pi P2 

where 

is the characteristic diffusion time, and 
1 



-TD) 



•>2] — 



\fWPlP~2 



P2Vl 
-PlP2 



(24) 



(25) 



(26) 



is the characteristic diffusion vector, which is orthogonal 
to |&). 

Now the problem is reduced to calculation of the ma- 
trix element (|25|), which contains the inverse matrix op- 
erator determined by integral collision operators. 



IV. EXACT AND LIMITING EXPRESSIONS FOR 
THE DIFFUSION TIME 

To derive an exact analytical expression for the un- 
known quantity ( j25| ) it is necessary to introduce a full sys- 
tem of orthonormal two-dimensional vectors \4>n) (here 
n = 1,2,3...) belong to the infinite-dimensional Hilbert 
space with a scalar product (18). The concrete choice 
of a system of basis vectors in many respects depends 
on convenience of calculations within the framework of a 
concrete physical problem (see, e.g. Refs. 5 and 8). In 
our problem it is convenient to take vector (|2^) as the 
first of them, and (|2^) as the second. The remaining vec- 
tors can be arbitrary (for example, such vectors can 
built on the basis of Sonin polynomials in classical case 
or Akhiezer-Aleksin-Khodusov polynomialsu in quantum 
systems), but should satisfy to the completeness and or- 
thogonality conditions 



^ \<l>m.)(<t>m\ = 1, {<t>m\<t>n) = &-n 



(27) 



In a constructed full system of vectors an exact expression 
for the diffusion time 



T D = -l[J+S 



(28) 



22 



can be reduced to the visual analytical form allowing sim- 
ple physical interpretation and providing a possibility to 
carry out the in-depth qualitative analysis of the obtained 
result in various limiting cases, corresponding to different 
mechanisms of equilibration in the system of quasiparti- 
cles: 



TD 



00 
n.m— 3 



(i+sy 



(29) 



Here the square matrices 2 and S contains the following 
matrix elements 



1511 



(<j>n\ J\<t>n 
(4>n\S\4> r , 



(30) 
(31) 



To obtain the formal result (|29|) I have used the relations 
© and @. 

The matrices ( |30|) and (|3l|) are infinite-dimensional 
and non-diagonal. Therefore the exact solution (|2^) does 
not allow to obtain the closed analytical expression for 
characteristic diffusion time. However it is necessary to 
emphasize, that the result ( p9|) contains explicitly not 
only the quantities responsible for interaction between 
quasiparticles of different types, but also the matrix el- 
ements appropriate to collisions between the quasiparti- 
cles within each subsystem. It allows to investigate vari- 
ous limiting cases, to find minimum and maximum values 
of diffusion time, to construct correct interpolation for- 
mulae and useful analytical models, to make calculations 
on computers for concrete physical systems. 

To carry out the detailed analysis of the formal solu- 
tion (^9|)' it is convenient to concretize slightly the vector 
basis \4> n ) for n > 2. Namely, let us choose the remaining 
vectors in the following form: 



1 



n<x+x = 



^ (Q) (Pi) 




\4>2a 



+2) 



(a) 




(32) 



(33) 



where a = 1, 2, 3, . . ., F^ a '(pk) is the system of properly 
chosen orthogonal polynomials with the norm Af% = 
(F(«)(p fe )| J F 1 ( Q )(p fc )) 1/2 , such that (F(«)(p fc )|p k ) - (it 
is clear that such polynomials cannot contain first degree 
of momentum). Such a choice (|32|), ( |33"| ) provide a pos- 
sibility to "separate" in some manner the component of 
the mixture. The matrix J + S takes now the following 
form 
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I v \ v \2 v \ 

(11) (11) (12) 



J + S 



y 21 



'21 



(21) (21) (22) 



V : 
o 



'\2 



(11) 



,(12) 



,(22) 



(34) 
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where I have introduced the following "higher" interac- 
tion rates: 



{a0) (F^( Pj )\C jk \FW(p k )) 



Jk 



and 



(af}) (F^( Pj )\C jk \FW( Pj )) 
for collisions of different quasiparticles, and 

3 rfritfP 



(35) 



(36) 



(37) 



corresponding to relaxation within each subsystem. The 
representation (^J) helps us to understand the explicit 
structure of the formal solution (^9|). This particularly 
provide a possibility to study the most important limiting 
cases and to reveal their physical matter. 

So, in case of infinitely fast establishment of an equi- 
librium between quasipaptieles of identical type (so-called 
complete control regimetSO) , when the strong inequali- 
ties 



(a/3) (a/3) (a/3) 



(j,k = l,2) 



(38) 



take a place, the second term in the brackets of general 
result (2^) vanishes and diffusion time ( pj| ) is given by 
the following simple formulae^ 



Jce) 
T D 



(min) 



where 



JO) 



-(Ckj 



1-22 



r (0)" 
12 



r (o)" 

'21 



(k,j = l,2;k^j). 



(39) 



(40) 



Here and thereafter the brackets (. . .) stand for a nor- 
malized average defined by the relation 

(A)k = ^-(p k \A\ Pk ), (k = 1,2). (41) 

According to the momentum conservation law C k j \p k ) — 
—Cj k \Pj) the relation between "basic" interaction rates 
in equation can be found as follows 



,(°)" 



P2 (0)- 
r 21 

Pi 



(42) 



Proceeding from hermiticity and negativity of opera- 
tors S and J with the help of well-known. Cauchy- 
Bunyakovsky inequality it is possible to showt3, that the 
following inequality is always valid 



\ (min) 



(43) 



where t^ 11 ^ is determined by eq. ( |39| ) 

In the opposite limiting case of extremely slow estab- 
lishment of equilibrium between identical quasiparticles, 
i.e. when 



r ) «^'Vjf ) , a* = 1,2), 



(44) 



the second matrix in cq. ( |34| ) can be neglected and the 
whole formula (|2^) can be converted so that the diffusion 
time td runs up to its maximum value 



^-)=-(^|J-i|0 2 ). 



(45) 



Comparing the main limiting results (|39|) and (|4g) we 
come to the important conclusion: The main qualitative 
difference between expressions for diffusion time in case 
of a fast and slow relaxation inside each subsystem of 
a mixture consists in a method of averaging of operator 
matrices corresponding to interaction between different 
subsystems. 

Very interesting situation can be realized when one, 
say first, subsystem equilibrates very slowly C\\ — ► 0, but 
relaxation in another one is extremely fast C22 — -> 00, so 
that 



,(«/3) 



(a/3) (a/3) 
l ^12 ' 



> V, 



(a/3) (a/3) 
i ^21 ' 



(46) 



In this case the results ( |29| ) and (34) can be converted as 
follows 



-(C12 )i - (C 2 i) 5 



(47) 



Let us pay attention to the principal difference between 
first and second terms in relation ( pl7| ) . Again, in the first 
term we average the time, while in the second one we av- 
erage the rate of interaction between quasiparticles of 
different species. According to the Cauchy-Bunyakovsky 
inequality the characteristic time (Cjk)^ 1 is always less 
than the time {C~ k l ) for any momentum dependence of 
Cjk- If the collision operator Cj k does not depend on 

the quasiparticles momenta, then 773 = Tq 11111 ^. If Cj k 
depends on momenta, then speed of equilibration of a 
system depends on how fast .the equilibration between 
identical quasiparticles iscJO. Such situation brightly 
manifests in a phonon-impuriton system of superfluid 
mixtures of helium isotopesca and in phonon systems in 
solidsE3. There is a two-stage mechanism of relaxation in 
these systems. At the first stage quasiparticles of the sec- 
ond type interact only with those quasiparticles of first 
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type, whose momenta correspond to maximum of colli- 
sion operator Ci 2 . At the second stage establishment 
of an equilibrium in a system is determined by the in- 
teraction of quasiparticles of 1st type with each other. 
This occurs because the quasiparticles of first type with 
minimum C12 prefer to interact with such quasiparticles 
within their subsystem, which are already at equilibrium 
with second subsystem. 

The significance of manner of averaging of collision 
operators can be illustrated on the simple example of 
lattice thermal conductivity in solids. The interaction 
rate between phonon and scattering center v p hi Cphi 
(say, point defect or impurity) is proportional to p^ h (so- 
called Rayleigh scattering). So, the integral (C~ hi ) p h di- 
verges at zero momentum and the corresponding relax- 
ation time tends to zero. This means that noncquilibrium 
long wavelength phonon simply " does not see" an impu- 
rity. At the same time the quantity (C p hi) P h is finite and 
leads to finite thermal current. So, the mechanism of 
equilibration of such system can be drawn as follows: at 
first all phonons come to quasi-equilibrium in their own 
subsystems, which corresponds to some stationary flux 
of phonons, then they begin to scatter on the impurities 
and point defects. As it can be seen from the results 
©, (|J), (H), (H|), and @, analogous competition 
mechanisms can occur in any two-component quasipar- 
ticle system. 

Calculation of (Cjk)j for particular physical systems 
does not meet any difficulties. However, to calculate 
(C~j})j we must inverse the collision operator Cjk, which 
is not a straightforward operation. As a rule, it can be 
done by replacing Cjk with some characteristic rate of in- 
teraction between quasiparticles i/jk = ^jfc(Pfc), and then 
by straightforward averaging the value VjkiPk)^ 1 , which 
is simply a multiplying operator. In the next Section I 
consider the most popular models in various problems 
associated with two-component classical gases or with 
condensed media whose transport properties are deter- 
mined by the processes in two-component quasiparticle 
systems. 



V. THE RAYLEIGH AND LORENTZ MODELS 
FOR TWO-COMPONENT GASEOUS MIXTURES 

Any theory claiming for a solution of some complicated 
special problems, should first of all be agreed necessar- 
ily with some fundamental results in the most simple 
limiting cases. In the kinetic theory of classical gases 
the diffusion in Lorentz gas (diffusion of light very rar- 
efied component in a gas of massive slow particles) and 
Rayleigh gas (diffusion of massive particles with small 
concentration in a light gas) traditionally are considered. 

Both classical models, i.e. the Rayleigh and Lorentz 
gases correspond in fact to the limiting case (^6|), but are 
more restricted in particles characteristics. Let us start 
with Lorentz model. This is the mixture of light compo- 



nent with very small concentration and the massive gas 
(buffer component). Let the light component be 1st one. 
Let me briefly review the classical approach to a problem 
of diffusion in such a mixturedo. The strong inequality 



ni/n 2 < 1 



(48) 



was required for light particles to interact only with a 
massive component, but not with one another. Large 
difference in masses of particles 



m 1 /m 2 -C 1 



(49) 



ensured an elasticity of scattering of particles of a light 
component on massive particles and large difference in 
thermal velocities of particles of different types. Thus, 
particles of the massive component can be treated as 
fixed and described by their equilibrium distribution 
function. With the purpose of calculation of a diffusion 
coefficient in such a system the strong inequalities ( f48| ) 
and ([l9|) are usually used for simplification of initial ki- 
netic equation, which can bfi—ceduced to the so-called 
Boltzmann-Lorentz equatiorotj. The consequent solu- 
tion of this equation gives explicit expression for a diffu- 
sion coefficient. 

Let us show, that the classical results for diffusion in 
Lorentz gas can be directly obtained from the general 
results given by eqs. (|2|) and (|34|). The operations below 
have certain value itself as a generalization of Rayleigh 
and Lorentz models on quantum gases of quasiparticles. 
Such generalization is not trivial because of impossibility 
to introduce the mass and conserved number density for 
some quasiparticles. So, we will formulate the problem 
in terms of normal densities. In this regard for Lorentz 
gas we replace two strong inequalities with the following 



P1/P2 < 1. 
Then we can rewrite vector |(/>2) 

1 




(50) 



(51) 



An equilibrium in massive component leads to the fol- 
lowing simple expression for operator matrix J given in 
general by eq. (|l6|): 



J = 



C12 




(52) 



Then, using the explicit expression ( p4| ) we can reduce 
the general result (|29| ) to the following form 



E {xi\Ci2\f^(pi)) (iK^+il^l^+i)!! -1 



a, 0=1 



a0 
-1 



(i^(pi)|C 12 |xi) • (53) 
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Vectors \<p2a+i) represent a complete set of orthonormal 
vector- functions in momentum space of first component. 
This allows us to rewrite relation (B3h as follows 



T D 



(oo) 
12 



(54) 



Further simplification can be achieved if we believe the 
massive particles to be fixed during a collision, so that 



vi - v 2 



vi 



(55) 



This particularly means that differential cross section de- 
pends on momenta of light particles only. Moreover, the 
momentum of a light particle can change only its direc- 
tion but not an absolute value. With these assumptions 
we can write collision operator Cyi as 



C12 = -\vi\a t ri2, 



(56) 



which is simply a multiplication operator and, therefore, 
it can be inversed without any difficulties. Here I intro- 
duced the transport cross sectionfl 



[1 — cos 9i2)da, 



(57) 



where 6\2 is the scattering angle, da is the ordinary dif- 
ferential cross section. As a result, using the relations 
(p5|)-(p7[) and defining_the diffusion coefficient in binary 
gas in usual mannema we come to the classical result 



Anni J a t 



(58) 



The obtained expression ( p8| ) coincides with analogous 
formulae obtained in Refs. 5 and 23 by simplification of 
the initial kinetic equation. 

Note that for deriving the formula (pq) in the frame 
of developed here approach it is enough to require ful- 
filment of strong inequality ( |50| ) for normal densities of 
components of a mixture, absence of equilibrium in first 
component, equilibrium in second component and elas- 
ticity of quasiparticlcs scattering. 

Now I shall consider one more classical— example, 
namely a diffusion in so-called Rayleigh gasE-3, that is 
diffusion of very rarefied massive gas in the light buffer 
component with large concentration. In other words, we, 
as well as in the previous case, have a mixture of light 
and massive components, but under opposite conditions. 
So we assume the concentration of a massive component 
to be so small that 



P2/P1 < 1. 



(59) 



Now I assume, that the light component of a mixture is 
already in equilibrium, while the particles of the massive 
component almost do not interact with one another. In 
this limiting case we again can use the result ( |53| ) with a 
replacement of subscripts 1 *-> 2. In this case, however, 



the relative velocity of particles from different compo- 
nents again is determined by a velocity of light particle. 
Therefore, the appropriate collision operator does not de- 
pend on momenta, so that in view of an orthogonality of 
vectors of selected basis its nongiagonal elements vanish, 



(x 2 |c 21 |F(«(p 2 ; 



0. 



03=1,2...). 



(60) 



Thus, proceeding from the relation (^3|) with a replace- 
ment of subscripts of components, we come to the rela- 
tion 



ART 



1 

3/0 2 



(p2|C 21 |p 2 )^r 2 ( 1 rl , 



(61) 



which determines diffusion time in a weak solution of 
massive component in an equilibrium light gas. Further, 
proceeding from a momentum conservation law, we ob- 
tain 



nin 2 f (o) 2 



(62) 



Expression (|62|) leads to the well known relation for a dif- 
fusion coefficient of massive particles in equilibrium light 
classical gasoa. Let me remark, that in Ref. 5 this re- 
sult was obtained by an indirect method with the use of 
Einstein relation between diffusion coefficient and mobil- 
ity. The use of the method developed here has allowed 
to obtain the formula (62) without engaging artificial ap- 
proaches, i.e. immediately from a general solution (|2S|), 
([m|). The developed approach allows to generalize the 
result (|6^) , as well as ( |58| ) , on quasiparticle systems with 
arbitrary dispersion law and statistics, in particular, on 
such systems, in which it is impossible to define notation 
of a mass in its classical sense. 



VI. KIHARA APPROXIMATION 

In the previous sections I considered only the most 
physically interesting limiting cases. However, in prac- 
tice the following problem can appear: How to calculate 
some dissipative coefficient more precisely in the interme- 
diate case, i.e. when the considered limiting situations 
do not take a place. Of course, the most straightfor- 
ward way is to compute it numerically using the formulae 
(|29|H|34|). In that case we are forced to restrict ourselves 
with some finite matrices in eq. (p9|). This is an ana- 
log of the Chapman-Cowling approximation in classical 
gaseous mixtures. In 1949 KiharacJ proposed another 
approximation, which is in fact simpler than Chapman- 
Cowling approximation. He simply proposed to neglect 
the nondiagonal integral brackets (see Ref. 2), which 
is exact for Maxwell molecules. Unfortunately, this ap- 
proximation is unproved until now and it can be partially 
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justified only by experience in numerical calculations for 
classical gases. In our theory such an approximation can 
be introduced by neglecting all the nondiagonal matrix 
elements in (|34| ) . After this procedure the inverse matrix 
(I + S) takes the form 

{i + sy 1 = 



ft™ 



t 



(11) 



t\ 





(22) 



V 



(63) 



where I introduced the diagonal characteristic times 

-l 



t 



(a/3) 



(a/3) . (a/3) 



(64) 



In view of eqs. (63) and (64) we can rewrite the result 
(29) in the following form 



D 



■(fa\c\<t> 



(fa\C\<f>s)(<h\C\<h 



(0 2 |C|^4){0 4 |C|^2) 

<0 4 |C|0 4 > 



(65) 



The series expansion fl65|) gives the correction to the lim- 
iting result (|39"|). In practice we should estimate numeri- 
cally an appropriate radius of convergence in series j65| ) 
and keep necessary number of terms. 

In general it is clear that the Kihara approximation 
works well when the momentum dependence of collision 
operators is weak enough. But sometimes the series like 
(65) have infinite radius of convergence. This indicates 
simply that it is better to start with the opposite limiting 
formula (45) as a zero order approximation. 



VII. CONCLUSIONS 

In the present work I present the general theory provid- 
ing a possibility to investigate diffusion processes in two- 
component gas of quasiparticles with arbitrary statistics 
and dispersion. The obtained main equations of the the- 
ory are correct for both systems with conserved and not 
conserved number of quasiparticles, which is mathemat- 
ically expressed in not equality or in equality of a chem- 
ical potential to zero, accordingly. The proposed theory 
can be generalized on the cases of classical and quantum 
gaseous mixtures with arbitrary number of components. 

To solve the formulated problem I start with the sys- 
tem of kinetic equations (Q) driving the evolution of 
corresponding distribution functions of components of a 
mixture. After a standard procedure of a linearization 
(0)' @ °f a kinetic problem, I have chosen the basis in 
infinite-dimensional two-parameter Hilbert space with a 



scalar product (Jig) selected so that operator matrix of 
collision integrals (|l5|)-(^7|) becomes Hermitian. As is 
known, the inverse matrix of collision integrals does not 
exist because of the moments of collision integrals. How- 
ever, by projecting on the nucleus of an integral operator 
of collisions it is possible to define somewhat inverse ma- 
trix. In the present work this procedure has been made 
with the help of introducing of projection operator ( [22] ) 
corresponding to conservation of total momentum of a 
quasiparticle system. As an outcome it allowed me to 
obtain the general expression ( p9| ) for characteristic time 
(see also eq. (pH)), determining speed of a diffusion re- 
laxation in a system. 

The obtained general result contains explicitly quanti- 
ties responsible for interactions between both quasipar- 
ticles from different components and between identical 
quasiparticles. It allows to analyze qualitative difference 
between mechanisms of equilibration of a whole system 
in various limiting cases. So, if the relaxation inside each 
component of a system is instantaneous, the diffusion 
time is determined by an inverse average of collision op- 
erator, describing interaction between quasiparticles of 
different types. In the opposite limiting case, when the 
equilibrium in a system occurs infinitely longly, the dif- 
fusion time is equal to an average of the inverse collision 
operator. The principal difference between these two lim- 
iting results can be easily understood using an example of 
phonon thermal conductivity of solids. The thermal con- 
ductivity in this case is simply a diffusion of phonons in a 
system of fixed "scatterers" (impurities, boundaries, de- 
fects etc.). If the phonons did not come yet to equilibrium 
with one another, the thermal conductivity is determined 
by average of inverse frequency of scattering of phonons 
on scatterers. In case of a long wavelength phonon, such 
frequency is proportional to the fourth degree of momen- 
tum (Rayleigh scattering). Therefore when averaging the 
integral simply diverges at zero momentum. That is the 
long wavelength phonon simply " does not feel" the scat- 
tered On the contrary, in case when the phonons came in 
an equilibrium with one another, the magnitude of diffu- 
sion time appears to be finite, because now the frequency 
is averaged, instead of time. This corresponds to the so- 
called two-stage mechanism of a relaxation. At first an 
equilibrium in phonon gas appropriate to some station- 
ary phonon flux is established. Then phonons scatter on 
scatterers. Just this process results in a finite heat flux. 

The general results, obtained in the work, give cor- 
rect expressions for diffusion coefficients of a light gas in 
a massive one (Lorentz gas) (pq) and a massive gas in 
a light one (Rayleigh gas) (j62|) . Earlier, these classical 
results were obtained only with use of artificial meth- 
ods based on simplification of the initial kinetic equation. 
Furthermore these two most popular model systems can 
be easily generalized to describe quasiparticle systems as 
well as classical gases. 
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